home *** CD-ROM | disk | FTP | other *** search
/ Graphics Plus / Graphics Plus.iso / general / fractal / kaos.lha / modellib / vanderpol_def.c < prev   
Encoding:
Text File  |  1990-01-24  |  983 b   |  58 lines
  1. int vanderpol_init()
  2. {
  3.     title_label = "Van der Pol's Eq";
  4.  
  5.     mapping_on = 0;
  6.     inverse_on = 1;
  7.     fderiv_on = 0;
  8.     enable_polar = 0;
  9.     enable_period = 0;
  10.  
  11.     var_dim = 2;
  12.     param_dim = 3;
  13.     func_dim = 2;
  14.  
  15.     (void) malloc_init();
  16.  
  17.     var_label[0] = "x";
  18.     var_label[1] = "y";
  19.     param_label[0] = "alpha";
  20.     param_label[1] = "beta";
  21.     param_label[2] = "omega";
  22.     func_label[0] = "t";
  23.     func_label[1] = "Undefined";
  24.  
  25.     param[0] = 1;
  26.     param[1] = 0;
  27.     param[2] = 1;
  28.  
  29.     var_i[0] = 0;
  30.     var_i[1] = 0;
  31.  
  32.     param_min[0]= -5; param_max[0]= 5;
  33.     param_min[1]= -5; param_max[1]= 5;
  34.     param_min[2]= 0; param_max[2]= 5;
  35.  
  36.     var_min[0]= -5; var_max[0]= 5;
  37.     var_min[1]= -5; var_max[1]= 5;
  38.  
  39.     f_p = vanderpol_f;
  40.     func_p = vanderpol_func;
  41. }
  42. /* Van der Pol's equation */
  43. int vanderpol_f(f,index,x,p,t,dim)
  44. int index,dim;
  45. double f[],x[],p[],t;
  46. {
  47.     double cos();
  48.     f[0] = x[1] - p[0]*(x[0]*x[0]*x[0]/3. - x[0]);
  49.     f[1] = -x[0] + p[1]*cos(p[2]*t);
  50. }
  51. int vanderpol_func(f,x,p,t,dim)
  52. double f[],x[],p[],t;
  53. int dim;
  54. {
  55.     f[0] = t;
  56.     f[1] = 0;
  57. }
  58.